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Abstract 

I review what we know about the donor stars in cataclysmic variables (CVs), 
focusing particularly on the close link between these binary components and the 
overall secular evolution of CVs. I begin with a brief overview of the "standard 
model" of CV evolution and explain why the key observables this model is designed 
to explain - the period gap and the period minimum - are intimately connected 
to the properties of the secondary stars in these systems. 

CV donors are expected to be slightly inflated relative to isolated, equal-mass 
main-sequence (MS) stars, and this "donor bloating" has now been confirmed 
observationally. The empirical donor mass-radius relationship also shows a dis- 
continuity at A/2 ~ 0.2Af© which neatly separates long- and short-period CVs. 
This is strong confirmation of the basic disrupted magnetic braking scenario for 
CV evolution. The empirical M2 — R2 relation can be combined with stellar mod- 
els to construct a complete, semi-empirical donor sequence for CVs. This sequence 
provides all physical and photometric properties of "normal" CV secondaries along 
the standard CV evolution track. 

The observed donor properties can also be used to reconstruct the complete evo- 
lution track followed by CVs, i.e. the mass-transfer rate and angular-momentum- 
loss rate as a function of orbital period. Such a reconstruction suggests that 
angular momentum loss rates below the period gap are too high to be driven 
solely by gravitational radiation. 

1. Introduction 

The last decade or so has seen tremendous progress in our understanding of 
cataclysmic variables (CVs), particularly as it relates to the evolution of these 
binary systems. Several of the key breakthroughs have been connected to the 
properties of the secondary stars in these systems. In this brief review, I will take 
stock of what we expect theoretically from and know observationally about CV 
donor stars. I will also explain how theory and observation can be powerfully 
combined to define a benchmark semi-empirical "CV donor sequence" and even 
to reconstruct the entire evolutionary path followed by CVs. 

2. The Evolution of Cataclysmic Variables: A Primer 

Figure [Li shows the orbital period distribution of CVs as of 2006. Two key fea- 
tures are immediately obvious: (i) the famous CV "period gap" between P gap .- — 
2 hrs and P ga p,+ — 3 hrs; (ii) the period minimum around P m in — 80 min. The 
most important requirement for any successful model of CV evolution is that it 
must provide a natural explanation for the origin and location of these features. 

In what has become the "standard model" of CV evolution, the period gap is 
interpreted as signalling a. switch in the dominant angular momentum loss (AML) 
mechanism. More specifically, the idea is that, above the period gap, CV evolution 
is driven mainly by "magnetic braking" (MB), i.e. by AML associated with a weak 
stellar wind from the donor star. This ionized wind is forced to co-rotate with 
the donor's magnetic field out to the Alfven radius, where the ram pressure in 



©2010 by Universal Academy Press, Inc. 



2 



t 1 1 — r 




Orbital Period (hrs) 

Fig. 1. Differential and cumulative orbital period distribution of CVs, based on data 
taken from Edition 7.6 of the Ritter & Kolb catalogue p] . Estimated values for the 
minimum period and the period gap edges are shown as vertical lines. The shaded 
regions around them indicate our estimate of the errors on these values. Figure 
reproduced from [2]. 



the outflow becomes equal to the pressure associated with the magnetic field. It 
turns out that the donor mass at the upper edge of the period gap corresponds 
roughly to the mass where the donor is expected to transition from a star with 
a radiative core to a fully convective object (Mi ~ 0.2 — 0.3 Mq, see Section 
3.1). The standard model thus posits that this transition effectively shuts down 
the magnetic field on the secondary and hence also disrupts MB. The physical 
justification for this idea is that the transition region between the radiative core 
and the convection zone - the so-called "tachocline" - is the location where the 
magnetic fields arc anchored in many magnetic dynamo models for low-mass stars. 

In the "vanilla-flavoured" standard model, MB ceases completely at the upper 
edge of the period gap, leaving only gravitational radiation (GR) to drive the 
further evolution of CVs. Why does such a switch in the AML rate produce a 
period gap? As it turns out, the answer is entirely associated with the properties 
of the donor star at the time of the switch and will be discussed in detail in Section 
3.3. For now, let us just note that the key difference between MB and GR is that, 
at least according to the standard model, the former produces much higher AML 
and mass-transfer rates, by a factor of > 10 in the vicinity of the period gap. 

Turning to the minimum period, this is often described as being associated with 
another transition of the donor, this time from a Hydrogen-burning star to a sub- 
stcllar object. The point here is that stars generally have a positive mass-radius 
index, whereas sub-stellar objects with masses below the hydrogen-burning limit 
(Mh — 0.07Mq) have a negative one. Thus, so long as M 2 > Mh, the donor 
radius, binary orbit, and orbital period are all expected to decrease in response 
to the mass loss the donor experiences, but all three quantities are expected to 
increase once Mi < Mh- We can therefore expect the condition Mi ~ Mh to set 
the minimum period a CV can reach. 
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One of the key points to take away from the discussion above is just how in- 
timately CV evolution is tied up with the properties of the donor stars in these 
systems. Consider: evolution above the gap is thought to be driven by a magnetic 
stellar wind from the secondary, the gap itself is thought to be associated with the 
transition of the secondary to a fully convective state, and the period minimum 
is set by the transition of the secondary into a sub-stellar object. In the following 
section, I will take a closer look at the physical link between CV secondaries and 
binary evolution. This will allow us to understand more precisely how the period 
gap is produced and also remind us that the period minimum is not necessarily 
associated with the stellar to sub-stellar transition of the donor. 

3. The Physics of CV Secondaries 

3.1. Fundamentals 

The radius of a Roche-lobc-filling star depends only the binary separation, a, 
and the mass ratio, q = M2/M1. A particularly convenient approximation for the 
Roche-lobe radius is [3] 
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which can be combined with Kepler's third law 
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to yield the well-known period- density relation for Roche-lobe- filling stars with 

If we are allowed to assume that donors are mostly low-mass, near-MS stars, we 
expect that their mass-radius relationship will be roughly 

R 2 /R Q = f(M 2 /M Q ) a (4) 

with / ~ a ~ 1. Combining this with the period-density relation immediately 
gives us approximate mass-period and radius-period relations for CV donors 

M 2 /M = M 2 /M Q ~ 0.1P orbM , (5) 

where P or b,hr is the orbital period in units of hours. This shows that, indeed, the 
period gap between 2 hrs and 3 hrs corresponds roughly to the expected transition 
of the secondary from partly radiative to fully convective, i.e. M 2 — 0.2 — 0.3M Q . 



3.2. Are CV Donors on the Main Sequence? 

So far, we have assumed that CV donors can be thought of as ordinary main- 
sequence stars. But is this actually true? 

The key question here is whether the mass loss the donor experiences should 
be expected to affect its overall stellar properties. The answer to this question 
depends on the competition between two time scales. First, there is the mass-loss 
time scale, 

r M2 - (6) 
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Fig. 2. The mass-radius relation for mass-losing stars for several ratios of Tth/r^, 
as indicated. Note that the notation in the figure legend is slightly different from 
that used in the text, but should be obvious: tkh = Tth an d tm = r A/ . Figure 
reproduced from [4]. 



which is the time scale on which the ongoing mass transfer reduces the donor 
mass. Second, we have the donor's thermal time scale, 

r th 10 8 (M 2 /A/ o )- 3 /Vs, (7) 

which is the time scale on which the donor can correct deviations from thermal 
equilibrium (TE). 

If mass loss is slow, in the sense that » Tth, the donor always has time 
to adjust itself to attain the appropriate TE structure for its current mass. It 
therefore remains on the MS, with a mass-radius index a ~ 1, and is essentially 
indistinguishable from an isolated MS star. 

Conversely, if mass loss is fast, i.e. t a/2 << Tth, the donor cannot adjust its 
structure quickly enough to remain in TE. Instead, the mass loss is effectively 
adiabatic. The response of low-mass stars with at least a substantial convective 
envelope to such mass loss is to expand, with a ~ —1/3. 

So which of these limits is appropriate for CVs? Neither, as it turns out. Let 
us take some typical parameters suggested by the standard model for CVs above 
and below the gap, say M 2 ~ 0.4 with M 2 ~ 1 X 1CT 9 and M 2 ~ 0.1 with 
M 2 — 3 x 10~ n , respectively. Plugging these values into Equations [5] and we 
find Tjj — Tth — 4 x 10 s yrs above the gap and r A ; /2 ~ Tth — 3 x 10 9 yrs below. 
Thus the thermal and mass-loss time scales are expected to be comparable for CV 
donors, both above and below the period gap. 

What does this mean for the response of the secondary to the mass loss it 
experiences? The answer is simply that the donor cannot shrink quite fast enough 
to keep up with the rate at which mass is removed from its surface. As a result, 
it is driven slightly out of thermal equilibrium and becomes somewhat oversized 
for its mass. This is nicely illustrated by Figure taken from [3], which shows 
how the donor mass-radius relationship depends on r t h I t m.^ , when this ratio is 
assumed to be constant along the evolution track. 
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3.3. The Importance of Being Slightly Disturbed 

So CV donors are almost, but not quite, in TE. Docs this slight deviation from 
TE actually matter? Yes it does. In fact, it is this deviation that is responsible 
for producing both the period gap and the period minimum. Let us see how this 
works. 

We will take as given, for the moment, that the period gap is "somehow" as- 
sociated with a sudden cessation of (or at least reduction in) MB a P or b — 3 hrs. 
But why should this produce a period gap in the CV population? 

Recall that the donor star is slightly out of thermal equilibrium - i.e. slightly 
bloated - as it encounters the upper edge of the period gap. Now, since mass 
transfer in CVs is driven entirely by AML, a sudden reduction in AML will also 
result in a sudden reduction in the mass-loss rate the donor experiences. This 
lower mass-loss rate cannot sustain the same degree of thermal disequilibrium 
and inflation in the secondary star. The donor therefore responds to this change 
by shrinking closer to its thermal equilibrium radius. However, this shrinkage 
almost immediately causes a total loss of contact between the stellar radius and 
the Roche lobe. The reason for this is that, in the semi-detached configuration, 
the stellar and Roche-lobe radii match extremely closely, to within |i?2 — Rl\ — H, 
where H/R2 ~ 1CP 4 is the exponential scale-height in the donor's envelope [5]. 
As a result, even a small reduction in radius (so long as it is greater than ~ H), 
will cause total loss of contact on a time-scale of roughly (H/R2)T t h ~ 10 4 yrs. 

The origin of the period gap in the standard model is now clear: a CV ap- 
proaches the upper edge of the gap with a slightly bloated donor star. The (as- 
sumed) cessation of MB at P or b — 3 hrs - associated, perhaps, with the transition 
of the donor to a fully convective state - then leads to a reduction in the mass-loss 
rate from the donor, which in turn causes the donor to shrink and lose contact 
with the Roche lobe altogether. The upper edge of the gap thus marks a cessa- 
tion of mass transfer in CVs. According to the standard model, CVs then evolve 
through the period gap as detached systems. During this detached phase, the 
binary orbit and Roche lobe continue to shrink, since there is still ongoing AML 
due to GR. However, provided the thermal relaxation of the donor is faster than 
the shrinkage of the Roche lobe, the donor manages to relax all the way back to its 
TE radius. In practice, this condition is met, so long as the AML rate is reduced 
by at least a factor 5-10 at the upper gap edge. The bottom edge of the period 
gap then corresponds to the location where the Roche lobe radius catches up once 
again to the TE radius of the donor. At this point, mass transfer restarts, and 
the system emerges from the gap as an active CV once more. 

How bloated must CV donors be to account for the observed size of the period 
gap? Since there is no mass transfer in the gap, the donor mass just above and 
below the gap must be the same, M2(P gap ,+ ) = M^Pgap,-)- From the period- 
density relation (Equation 13]), we then get 

R2{Pgap,+ ) 
R2(Pgap,-) 

We also know that the donor at the bottom edge is in or near equilibrium, so 
we conclude that donors at the upper edge of the period gap must be oversized by 
~ 30% relative to equal-mass, isolated MS stars. 

Let us now take a closer look at the minimum period for CVs. As it turns out, 
this, too, is closely connected to the properties of the donor stars in these sytems. 
If we combine the period-density relation (Equation [3|) with the simple power-law 
approximation to the donor mass-radius relation (Equation [4]), we find 



p 



p 



gap. 



1.3. 



(8) 



P~ 2 rv A/f 1_3a 
r orb KM 2 



(9) 



6 




M a (M ) 



Fig. 3. The mass-radius relation of CV donor stars, based on the data presented in 
[6]. Superhumpers are shown in black, eclipsers in red (this may not be visible in 
the printed version of this manuscript; see Figure 5 for a clearer view of the same 
data in this case). Filled squares (circles) correspond to short-period (long-period) 
CVs, crosses to likely period bouncers. The parallelograms illustrate typical errors. 
Open symbols correspond to systems in the period gap or likely evolved systems. 
The solid lines show the optimal broken power-law fit to the data. The dotted line 
is a theoretical mass-radius relation for MS stars [7J. Figure adapted from [2J. 



Differentiating this logarithmically yields a simple expression for the orbital period 
derivative, i.e. 

P orb _ 3a - 1 M 2 , , 

Porb 2 M 2 ■ [ > 

Since the period minimum must correspond to P or (, = 0, Equation [TU] tells us that 
Pmin occurs when the donor has been driven so far out of thermal equilibrium 
that its mass-radius index along the evolution track has been reduced from its 
near-MS value of a ~ 1 to a = 1/3. So, as already noted above, P m i n does 
not necessarily have to coincide with the orbital per iod a t which the donor mass 
reaches Mjj. In fact, recall that we noted in Section that, for any donor with 
at least a substantial convective envelope, the mass-radius index in the limit of 
fast (adiabatic) mass-transfer is a ~ —1/3. Thus the period evolution of a CV 
can in principle be made to turn around at any donor mass, provided only that 
mass loss becomes sufficiently fast compared to the donor's thermal time scale. 
The significance of Mh in this context is that period bounce becomes inevitable 
when the donor reaches this limit. This is because sub-stellar objects are out of 
TE by definition and respond even to slow mass loss by increasing in radius, i.e. 
a < 0. And, in practice, P m in does, in fact, correspond roughly to M 2 — Mh- 

4. Are CV Donors Observationally Distinguishable from MS Stars? 

We have seen in the previous section that we expect CV donors to be significantly 
larger than equal-mass, isolate MS stars. Can we actually observe this difference? 

Yes, we can. Figure [IT] shows the empirical mass-radius relationship of CV 
donors, as first constructed by [6]. The figure shown is actually taken from [2], 
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but this is only a minor modification of the original relation in [5] and is still based 
on the same data. 

Figure [XI represents a major break-through in our understanding of CV donors 
and evolution. First, it definitively confirms the theoretical expectation that CV 
donors are larger than ordinary MS stars, both above and below the period gap. 
Second, the donor mass-radius indices above and below the gap are just what the 
doctor ordered: they are less than the MS-based (TE) values, but greater than 
the critical value of a = 1/3. Third, the period bouncer regime is very poorly 
constrained by the data, but we do find M2 ~ Mh near P m in, as well as a < 1/3 
for the lowest mass donors below this limit. Again, this is nicely in line with our 
theoretical expectations. 

However, the single most important aspect of Figure [23 is that it reveals a 
discontinuity in the donor radii at M% ~ 0.2Mq that neatly separates long-period 
CVs from short-period ones. Moreover, on the low-mass (short-period) side of 
the discontinuity, i? 2 is very close to the MS (thermal equilibrium) radius for this 
mass. This is by far the strongest evidence to date that the basic "disrupted MB'' 
scenario for CV evolution is fundamentally correct. 

5. A Complete Semi-Empirical Donor Sequence for CVs 

Virtually all observable properties of CV donors depend on just three physical 
parameters: M2, R2 and T e ff^- For example, the total luminosity of the donor 
is Li = AwR20~Tgfj! 2 , its spectral energy distribution depends primarily on T e //,2 

and log 52 = log (GM2/R2), and, with the SED fixed, its flux in any particular 
wave-band depends only on R2. Now the empirical M2 — R2 relation in Figure [3~1 
gives us a unique relationship between two of these three key donor parameters. If 
we could find just one additional relationship between either of these two param- 
eters and T e ff^2, we would essentially know all there is to know about "ordinary" 
CV donors. 

As it turns out, there is such a relationship, albeit a theoretical one. The key 
insight is that the low-mass donors we care about all have large convective en- 
velopes. As a result, their effective temperature is almost completely independent 
of luminosity and instead depends only on mass [51 12] ■ Thus CV donors are ex- 
pected to obey a standard MS M 2 — T e ff 2 relationship, which can be taken from 
state-of-the-art stellar models. By combining the empirical M2 — R2 relation with 
the theoretical M2 — T e fj.2 one, we can then construct a complete, semi-empirical 
donor sequence for CVs that defines all physical and photometric properties of 
CV secondaries as a function of P or b- 

This program was carried out in [2\ , and some key aspects of the resulting donor 
sequence are shown in Figure [4] This sequence is useful for at least two reasons: 
first, it helps to define precisely what we mean by a "normal" CV, i.e. it provides 
a useful benchmark for theoretical and observational studies of CV evolution. 
Second, it provides an immediate estimate of the expected donor brightness for 
any CV in any desired waveband, given only an estimate of the system's P or b- 
This can be used, for example, to set limits on distances to CVs, via the method 
of photometric parallax. 

6. Reconstructing CV Evolution from Donor Properties 

The ultimate goal of essentially all work on CV evolution is to determine the 
AML rate as a function of P or b; i.e. to find the correct form of J(P or b), or, 

equivalently, A'hiPorb)- The reason this is the "holy grail" is that such a recipe 
would allow us to calculate/predict almost everything there is to know about 
CVs, from the complete set of binary properties of individual systems at given 
Porb, to the population properties of large samples of CVs (e.g. their orbital 
period distribution). In reality, there would, of course, still be complications - 
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Fig. 4. Physical and photometric parameters along the semi-empirical CV donor se- 
quence. The left-hand column shows physical donor parameters as a function of P or b, 
while the middle and right-hand columns show optical and near-infrared absolute 
magnitudes, respectively. The dashed lines show the expected absolute magnitudes 
of the accretion heated white dwarf. Figure reproduced from [2]. 



for example, we would still have to worry about selection effects [TB] - but it is 
nevertheless true that "understanding CV evolution" is broadly synonymous with 

"knowing j{P or b) and/or M 2 (P orb )" . 

It may seem that we have already extracted a lot of information from the 
empirical M 2 — R2 relation in constructing our semi-empirical donor sequence. 
However, we can actually push things even further and use the M 2 — i?2 relation 
to reconstruct the full evolution path of CVs, i.e. to determine M2(Porb) and 

hence J(P or b). It is actually easy to see that this should be possible: after all, the 
degree of thermal disequilibrium and radius inflation a donor experiences increases 
with increasing mass- loss rate (see, for example, Figure [2]). Thus R2{P rb) is a 
direct tracer of M 2 (P or b)- 

The great advantage of donor-based M 2 {P or b) and J{P or b) determinations is 
that they are very likely to trace the true secular (i.e. long-term) rates. This is 
primarily because the time scale on which the donor can adjust its radius is long 
compared to th e averaging time scales inherent in essentially all other methods 
[TU1 [TT1 Ir21 [T5] . The main difficulty, on the other hand, is that one has to care- 
fully correct for other effects that might cause (or masquerade as) donor-bloating, 
such as stellar activity, imperfect stellar models, tidal/rotational deformation and 
irradiation. 

Figure ED (taken from [14]) shows the results of an attempt to deal with these 
complications and derive a complete, purely donor-based evolution track for CVs. 
The observed M2 — R2 points are the same as in Figure EU but the figure now 
also shows two self-consistent CV evolution tracks superposed on this data set. 
The thin black line shows the track predicted by the standard model, in which 
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Fig. 5. The observed donor M2 — R2 relation compared to the theoretical relations 
predicted by two self-consistent CV evolution tracks. The same symbols are used 
for the data points as in Figure 3, except that eclipsers are now indicated by open 
diamonds. The black dash-dotted line shows the M2 — R2 relation predicted by the 
standard model for CV. In this model, AML above the period gap is assumed to 
follow a standard "RVJ-like" MB law with 7 = 3 [15], while AML below the gap is 
assumed to be driven solely by GR. By contrast, the red line shows the best-fitting 
M2 — R2 relation, if we allow the strength of AML to deviate from the standard 
prescription. Figure adapted from |14j . 



MB is assumed to be at full strength above the gap, and switches off completely 
at P gaPl +- By contrast, the thick red line shows the optimal fit to the data, which 
requires approximately 0.65 x the standard MB rate above the gap, but 2.5 x the 
GR-drivcn AML rate below. 

Depending on one's point of view, this best-fit model is cither a minor modifica- 
tion of the standard model (in the sense that is still a straightforward implementa- 
tion of the basic disrupted MB scenario) or a major departure from conventional 
CV wisdom (in the sense that it requires AML rates significantly in excess of GR 
below the period gap). What is clear, however, is that the revised model fits the 
donor data substantially better than the standard one. 

The revised model may also resolve two long-standing problems in our under- 
standing of CV evolution. First, the minimum peri od predicted by the stand ard 
model is considerably shorter than is observed [17]. By contrast, we find in [14] 
that the revised model does an excellent job of matching the observed location 
of Pmin [II]- Second, observations suggest that the number of long-period CVs 
relative to short-period, pre-bo unce CVs i s higher than predicted by the standard 
model, by at least a factor of 3 [T6] [T9l [20] . In the revised model, the combination 
of enhanced AML rates below the gap and reduced rates above the gap increases 
the predicted ratio by just this factor [14]. 
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7. Summary and Conclusions 

Our understanding of CV secondary stars, as well as their relation to CV evolu- 
tion, has improved dramatically over the last decade or so. Perhaps most impor- 
tantly, we now know empirically that CV donors are oversized relative to equal- 
mass MS stars, and also that their mass-radius relation has a discontinuity at 
Mi ~ 0.2M© that separates short-period and long-period CVs. All of this is 
strong confirmation of the basic disrupted MB scenario for CV evolution. 

By combining the observed M-x — R2 relation with a theoretical M2 — T e f f.2 one, 
we have also been able to construct a complete, semi-empirical "donor sequence" 
for CVs that provides all physical and photometric parameters of CV secondaries 
as a function of only P or b- 

Finally, we can even use the observed M 2 — R2 relation to reconstruct the entire 

evolutionary path of CVs, that is to say J(P or b) and M2(P or b)- This is possible 
because the degree of donor inflation is a direct measure of its mass-loss rate (and 
hence of the AML rate from the system). A first attempt to implement this idea 
suggests that the MB-driven AML rate above the period gap is slightly lower than 
is usually assumed, but that an AML rate in excess of GR (by a factor of ~ 2.5) 
is required to match the observed mass-radius relation. This revised model may 
also resolve two long-standing problems in CV evolution: the mis-match between 
the observed and predicted location of P m in, and the higher-than-expected ratio 
of long-period CVs to short-period, pre-bounce CVs. 
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